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We present a method for the experimental measurement of nonclassicality witnesses and demonstrate its 
application. Our proposal only requires the coherent displacement of the initial state, which can be achieved by 
overlapping the latter with a coherent state at a beam splitter, and subsequent photon-number resolved detection. 
This setup allows a complete test of nonclassicality of an arbitrary quantum state. The role of the quantum 
efficiency as well as statistical and systematic uncertainties are discussed. Finally, the scheme is demonstrated 
for a realistic example. 
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I. INTRODUCTION 

Nonclassicality describes the difference between classical 
and quantum physics. For the harmonic oscillator, its formal 
definition is based on the coherent states |a), which are seen 
as the closest analogues to the classical oscillation |Q]. Re- 
markably, any quantum state p can be formally written as a 
statistical mixture of (classical) coherent states ||2[3], 

p = J d 2 aP(a)\a){a\. (1) 

The function P(a) is the so-called Glauber Sudarshan P func- 
tion, playing the role of the probability density on the set of 
coherent states. However, it may violate the requirements for 
a probability density, such that it has to be considered as a 
quasiprobability. A state p is referred to as classical, if the 
P function still fulfills all properties of a probability density. 
Otherwise, the state is called nonclassical |4|. 

However, this definition cannot be applied in most practi- 
cal situations, since the P function can only be understood 
as a generalized function, having singularities of (5-type and 
worse. Therefore, several derived criteria for nonclassicality 
have been developed. On the one hand, one may look at differ- 
ent quasiprobabilities. For instance, negativities of the Wigner 
function are indicators for nonclassicality of a state |5] |6), 
but they do not appear for all nonclassical quantum states. 
The concept of nonclassicality filters and quasiprobabilities 
is more powerful, since it provides a simple and complete 
characterization of nonclassical effects Q. First experimental 
demonstrations have been presented in Refs. ll8j {T0l . 

On the other hand, one can examine the expectation value 
of a so-called nonclassicality witness W for a given quantum 
state lfm[T2ll . If this expectation value is not compatible with 
expectation values of classical states, then the state must be 
nonclassical. This criterion has been used to formulate condi- 
tions for matrices of moments [ 13 1 or conditions for outcome 
probabilities lfl4ll . 

Recently, the relation between nonclassicality witnesses 
and quasiprobabilities has been elaborated [15|, providing 
a unified view on all kinds of nonclassicality tests. More- 
over, it has been shown that sets of nonclassicality witnesses, 
which enable one to detect nonclassicality of an arbitrary 
quantum state, can be parameterized by three real numbers 
only. Explicit examples for such sets of witnesses have al- 



ready been given. In the present work, we show in detail how 
one can measure the expectation values of these witnesses, 
by using coherent displacement of a quantum state and sub- 
sequent photon-number resolved detection. The experimen- 
tal setup is similar to the unbalanced homodyne measurement 
scheme 1 16 1, but adapted for our particular application. This 
enables us to determine the expectation values at the quantum- 
mechanical level of uncertainty ifTTl . which can be used for 
the significant verification of nonclassical effects. 

The paper is structured as follows. In Sec. II, we briefly re- 
view the construction of complete sets of nonclassicality wit- 
nesses. Then, we give the theoretical description of the dis- 
placement and measurement in Sec. Ill, and discuss the role 
of statistical and systematic uncertainties. Sec. IV is dedicated 
to an illustrative example. A brief summary and some conclu- 
sions are given in Sec. V. 

II. THEORETICAL BACKGROUND 

A nonclassicality witness is an operator W, whose expec- 
tation value is nonnegative for all classical states, 

(a\W\a) > 0. (2) 

Clearly, if one examines a state p for which the expectation 
value of W is negative, then the state p is nonclassical: 

Tr-j^VK} < => p nonclassical. (3) 

It has been shown in [15| that for any nonclassical state p, 
one can verify the nonclassical character by taking a set of 
nonclassicality witnesses of the form 

W w (a) = D(a)W w D(-a), (4) 

where D(a) = e aat ~ Q a is the coherent displacement opera- 
tor, depending on the complex amplitude a, and the operator 
W w is defined as 

W w — w 2 :uj*(wa , wa)us{wa\ wa) : . (5) 

Here : • : denotes the normal ordering procedure, and w is a 
real, nonnegative width parameter. 

The function ui(a* , a) can be chosen fixed. It only has to 
satisfy some weak properties, namely the Fourier transform of 
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\w{a* 



a 



shall be a nonclassicality filter fulfilling the con- 



ditions as introduced in |7|. Roughly speaking, this means 
that u(a*,a) is a real continuous function, and the Fourier 
transform of \co(a* ,a)\ 2 is decaying faster than any Gaussian 
function. A typical example of a witness operator was derived 
in 1151. 



W w (a) 



[Ji(wy/(at -q*)(a-a))] s 
4(fit — a*) (a — a) 



(6) 



where J\(x) is the Bessel function of first order. Hence, a 
universal witness applying to any quantum state can even be 
given in a closed analytical form. The price we have to pay 
is that such a witness is based on a nonclassicality filter with 
compact support. This means that in general it does not test 
the complete properties of the quantum state under study for 
a chosen w parameter. However, by increasing the value of w, 
that test converges to a complete one. 

To overcome this problem, we can choose a witness based 
on autocorrelation filters 0Q21- It tests the full state for any 
choice of w, but can be represented only via a numerical inte- 
gration. This, however, is not a serious problem since one can 
store the result of the numerical integration to further apply it 
for the direct handling of the experimental data. 



III. MEASUREMENT OF NONCLASSICALITY 
WITNESSES 

A. Experimental setup 

In the following, we restrict ourselves to a function of the 
form u){a' , a) = u>(h), where n — a) a is the photon number 
operator. In this case, the set of witnesses Q can be written 
in the form 



^(q) = w 2 :w\w 2 n{a))w{w 2 h{a)) :, 



(7) 



with fi(a) = D(a)hD(—a) being the displaced photon num- 
ber operator. In [16|, a procedure has been proposed to mea- 
sure the statistics of such operators. 
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Figure 1. Experimental scheme for the detection of a nonclassicality 
witness. The beam splitter displaces the input state p, the photodiode 
PD records the photon number statistics with quantum efficiency r\o- 
The expectation value of the witness W u (a) can be calculated from 
Eq. (TTTJ. 

Let us look at the sketch in Fig. [T] The state of interest 
is overlapped at a beamsplitter with a coherent state of am- 
plitude ta/r, where t, r are the transmissivity and reflectivity 



of the beamsplitter. The output state is measured by a pho- 
ton number resolving detector with quantum efficiency -qu . It 
has been shown that the recorded statistics corresponds to the 
statistics of the operator h(a) of the input state p, measured 
with an effective quantum efficiency r\ = t 2 r\D- Therefore, 
with this setup, we obtain probabilities p n (a, if) for observing 
the n-th outcome of the measurement of the photon number 
operator h(a). 

Now, we can use the witness operator W w in its spectral 
decomposition into photon number states |n), 



n=0 



with the matrix elements 



Qw,n — (n\ '■ w\w 2 h)uj(w 2 n) : \n). 



(8) 



(9) 



Note that these numbers are independent of the coherent dis- 
placement a. All off-diagonal matrix elements are equal to 
zero, since we assumed that W w does only depend on the pho- 
ton number operator. Analogously, we can find the spectral 
decomposition of the displaced witness as 



W w (a) 



OO 

E 

n=0 



n w ,n\n{a)}{n{a)\, 



(10) 



where |n(a)) = D(a)\n) denote the eigenstates of the dis- 
placed photon number operator. Now let us calculate the ex- 
pectation value of the operator in Eq. ( fT0| ) with respect to an 
input state p v , which suffered detection losses i]. The quantity 
is given by 



(W w (a))v - Tr{p v W w (a)} 



OO 

E 

n=0 



(11) 



where p n (a, rj) — (n(a)\p v \n(a)) is exactly the statistics 
of the photon number which is measured in the described 
scheme. Together with the theoretically known coefficients 
f2 w , n , this equation allows us to estimate the expectation value 
of the witness. 

Equation (JTTJ is a central result of this manuscript, which 
enables us to look for negative expectation values of any wit- 
ness W w (a). We emphasize that for verification of nonclassi- 
cality, it is sufficient to observe negativity of the expectation 
value for a single complex amplitude a. In this sense, non- 
classicality can be tested locally in phase space. Conversely, 
to demonstrate classicality of a state, one has to examine the 
witnesses of all complex amplitudes a. 

A strong point of our method consist in the following. In 
typical experimental scenarios one has at least a clear conjec- 
ture about the state one is going to produce. Hence on can 
easily chose the amplitude a to be optimal for the state one 
expects to be created. When this strategy is successful, the 
test only requires the optimization of one real parameter w. A 
failure of the test proves the failure of the conjecture. In such 
cases one has still the option to vary the amplitude a. In any 
case, the test procedure only requires photon-number resolved 
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measurements in order to estimate the witness for an arbitrary 
amplitude a. We note that our witness test can be replaced 
by any other nonclassicality criterion, which is necessary and 
sufficient for photon number distributions. However, they can 
be much more elaborate than the calculation of a simple ex- 
pectation value |T8l . 

One can minimize the losses introduced by the beamsplitter 
by choosing a transmissivity close to one. However, theoreti- 
cally the quantum efficiency is not important for the examina- 
tion of nonclassicality, since it does not destroy nonclassical 
effects, but may only decrease the significance of their detec- 
tion. More precisely, due to the equivalence of nonclassicality 
witnesses and quasiprobabilities 1031 . we can adapt the con- 
siderations in [8 1 and find that the quantum efficiency only 
leads to rescaled parameters w, a of the witness, 



r,Tr{p v W w (a)} = Tv{pW^ w (a/ ^j)}. 



(12) 



Since it is known that the nonclassical effects become more 
pronounced with increasing w, a lower quantum efficiency 
may make their significant verification more difficult. How- 
ever, as a matter of principle, it does not change the property 
of a quantum state of being nonclassical or not. Therefore, 
for any positive quantum efficiency it is possible to verify 
nonclassicality if the number of recorded data points is suf- 
ficiently large. 



B. Statistical uncertainties 

Since the photon number probabilities p n (a,rj) are mea- 
sured experimentally, they are subject to a statistical uncer- 
tainty, stemming from the finite number of measurements. 
Typically, they are calculated from the number of occurrences 
N n of a specific photon number n, divided by the total number 
of measurements N. Since the N n follow a multinomial dis- 
tribution, with probabilities p n {a, ry) for each outcome, their 
covariance matrix is given by [ 19 1 



C(N m ,N n ) = Np n (a,T])(S mn - p m {a:, t)). 



(13) 



Consequently, the covariance of the estimated probabilities 
equals to C(N m , N n )/N 2 . Now, the expectation value (11 
is a linear function of the photon number probabilities. There 
fore, its variance can be exactly estimated as 



1 



OO 



Vta((W w (a))„) = ^ C{N m , N n )fl w , m ft w>n . (14) 

m,n— 



Obviously, the evaluation of statistical errors is straightfor- 
ward and very simple. Moreover, we note that the nonclas- 
sicality test is performed on the quantum-mechanical level of 
uncertainty, since only measurements of a single observable 
n(a) are required ifTTl . 



C. Systematic error estimation 

In practice, we cannot measure an infinite number of pho- 
ton number probabilities p n . However, we also know that all 



states with a finite maximum photon number are nonclassical. 
Therefore, it is necessary to estimate the error which arises 
due to the restriction of the state in Fock basis, in order to 
evaluate if a given state is nonclassical or the observed nega- 
tivities are only due to the truncation of the Fock basis. 

Let us assume that we can measure the probabilities of the 
Fock states |0) to \N). Therefore, we estimate the expectation 
value of the witness by 



N 



p n (a,r))S} W! ,, 



(15) 



The corresponding error therefore is given by 



AW w (a,p) = Pn(a,1j)^w,n (16) 

n=N+l 

For a theoretically known state p, it is not difficult to evalu- 
ate this quantity. However, we do not want to use assump- 
tions about the state we examine. Therefore, we have to find 
a suitable alternative. Here, we ask for the maximum of the 
deviation AW w (a, p) for all classical states, 



AW W , 



max[AW w (a, p c i,)] 

Pel. 



(17) 



If we can estimate this number, we know that if the state is 
classical, its expectation value of the nonclassicality witness 
is bounded from above by 



(W w (a)) ieah < (W w (a)) 

I 



AW W 



cl.- 



(18) 



If the right side of this equation still shows significant negativ- 
ities, the state cannot be classical, but must be a nonclassical 
one. 

As one can imagine, the calculation of the systematic er- 
ror is more involved, such that we give the details in the Ap- 
pendix. For practical reasons, we do not estimate the maxi- 
mum with respect to all classical states, but restrict ourselves 
to states for which the probability for counting more than N 
photons is less than some threshold 1 — p r , with p T being the 
total probability of all resolved photon numbers. This restric- 
tion is based on the assumption that the recorded data already 
gives sufficient information on the photon statistics, and the 
photon counts which cannot be resolved by the detector are 
negligible. It turns out that the systematic error is independent 
of the point a, therefore we already omitted the argument in 
AW w , c i. . Finally, we show that the maximum can be found as 
the extremal point of a function of two real parameters, repre- 
senting amplitudes of coherent states, cf. Eqs. ( f52| ), ( |53| l. This 
makes the calculation feasible for practical applications. 



IV. EXAMPLE: SINGLE PHOTON ADDED THERMAL 
STATE 

Let us consider a single photon added thermal state 
(SPATS) ll20l with a mean photon number n = 0.8, and a 
preparation efficiency 77 prop . = 0.5. It is known that this state 
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exhibits neither squeezing nor sub-Poissonian statistics, aris- 
ing from Eq. <J3J) with W =:(£ — (x)) 2 : and W =:(n — (n)) 2 - 
respectively iflll . Therefore, these two well-known nonclas- 
sicality witnesses are not able to verify nonclassicality of this 
particular state. It is also important that under the chosen con- 
ditions the Wigner function is positive semidefinite ETl . and 
the characteristic function does not show nonclassicality of 
first order Ell . Therefore, it is a nontrivial task to witness the 
nonclassicality of the state under study. 

We assume that the quantum efficiency of the detector 
equals to t)d — 0.5, leading to an overall efficiency of rj = 
0.25. Moreover, let us suppose that the detector may resolve 
the photon number statistics for up to TV = 15 photons. The 
initial beam splitter has a transmissivity t = 0.99. We fur- 
ther assume that we perform 10 5 measurements per coherent 
displacement a. The probability of finding not more than 15 
photons equals to p r = 0.995 for coherent displacements with 
| a | < 3. As a nonclassicality witness, we choose the example 
given in Eq. |6]). Then, the matrix elements |9]l are given by 

= w 2 " (_ w 2 /4) m / 2to + 2 \ ml 

W " n 16 i( m + 1 V-] 2 V m ){n-m)V 

Since SPATS are phase symmetric and similar to the photon, 
we may test the witness at the origin of phase space, i.e. for 
a = 0. 
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Figure 2. Dependence of the expectation value of the witness W„, (0) 
on the filter width, for a SPATS with mean thermal photon number 
n — 0.8 and overall efficiency r\ — 0.25. The blue shaded area 
corresponds to a statistical error of one standard deviation, the red 
shaded area shows the systematic error due to the truncation of the 
Fock space. 

In Fig. [2] the expectation value of Ww(0) is shown in de- 
pendence of w together with the statistical and systematic un- 
certainty. We note that the systematic error strongly depends 
on the value of w. At some points between w = 4 and w = 5, 
it plays an important role, while it can be neglected for other 
values of w. For sufficiently large w, the expectation value be- 
comes clearly negative. The optimal significance is achieved 
for w = 4.2. For this choice of the filter widt h our result 

is (W w (0)) + AW w , ct = -6.lJvar((W w (0))), which is a 



significance of more than six standard deviations. 
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Figure 3. Dependence of the expectation value of the witness on 
the coherent displacement a, for a SPATS with mean thermal photon 
number n = 0.8 and overall efficiency rj = 0.25, and an optimal 
filter width w = 4.2. The blue shaded area corresponds to a statistical 
error of one standard deviation, the systematic error is negligible for 
this width parameter. 

The dependence of the expectation values of W„(q) on the 
coherent displacement a is shown in Fig. [3] The width w is 
chosen to be the value for optimal significance. The system- 
atic error is negligible in this case. We observe the largest 
negativity for a = 0. There, nonclassicality is verified with 
sufficiently high significance. We note that the curve also 
shows the course of the nonclassicality filtered P function of 
the SPATS, cf. EJUS] 

Let us conclude this section with a brief discussion of the 
photon number resolution which plays a key role for the 
method under study. Whenever for the light source under 
study a properly mode-matched local oscillator (LO) is avail- 
able, the photon number resolution can easily realized by bal- 
anced homodyning with a phase randomized LO ll23l . For the 
considered example of a SPATS this method could be imple- 
mented EH, so that the required photon number resolution of 
up to 15 photons is easily realized. For other radiation sources 
a proper LO may not be available. For such situations nowa- 
days different realizations of photon number resolving detec- 
tors are under study 1 24-29 1 . The theory of such detectors 
shows that the measured statistics only slowly converges to 
that of a truly photon number resolving detector [30]. How- 
ever, one may directly use a measured sub-binomial statistics 
to certify nonclassicality [31|. The further development of 
such direct methods for the general witnessing of nonclassi- 
cality is a demanding task for future research. 

V. SUMMARY AND CONCLUSIONS 

We have analyzed a simple experimental setup to measure 
the expectation value of an arbitrary nonclassicality witness. 
The coherent displacement a can be chosen by the experi- 
mental setup, only the width parameter w has to be introduced 
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by a suitable post-processing procedure. Its optimal value de- 
pends on the quantum state under study and on the amount of 
available recorded data. It can be properly adjusted during the 
rather elementary computational data analysis. 

The applicability of our method has been simulated under 
realistic conditions for a single-photon added thermal state. 
For the analyzed quantum state, the commonly used nonclas- 
sicality criteria fail to identify its nonclassicality. The state is 
neither sub-Poissonian nor quadrature squeezed, the Wigner 
function is positive semidefinite, and the characteristic func- 
tion does not display nonclassicality directly by exceeding the 
boundaries set by its vacuum noise level. Even under such 
complicated conditions and after taking systematic and statis- 
tical uncertainties into account, our witness approach clearly 
verifies the nonclassicality with a good significance. 
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APPENDIX: SYSTEMATIC ERROR ESTIMATION 

Here, we show how to calculate the systematic error 



AW W h = max 



,n=N+l 



(20) 



All classical states can be commonly described by a nonneg- 
ative P function, therefore we take the optimization with re- 
spect to all P functions admitting an interpretation as a classi- 
cal probability density. Clearly, the objective function is linear 
in P(a). Therefore, we have to solve a generalized linear op- 
timization problem, with constraints 



P{a) > 0, J P(a)d 2 a = 1. 



(21) 



We solve this task in several steps. First, we formulate 
the optimization problem and transfer it into a more simple 
form. Second, we discretize the problem and apply standard 
knowledge about linear optimization problems to uncover the 
structure of the solution. Last, we argue why this structure 
also holds in the continuous case, and present a simple equiv- 
alent constrained optimization problem, which only requires 
to maximize with respect to two real variables. 



1. Optimization problem 

Let us examine the objective function for coherent states 
p c i. = 



n=N+l 



Here we already inserted the photon number distribution 
p n (a, v)- Obviously, the objective function only depends on 
the value of \a — (3\^/fj. Therefore, we may write 

AW w (a,\l3},ri) = AW w (\a-p\^Tj,\0),l). (23) 

For arbitrary states, characterized by a P function P(f3), the 
systematic error can be obtained via 

AW w (a,p,r))= [ d 2 pP(/3)AW w (a,\P),r,) 

d 2 pP(p)AW w (\a- p\^/rj,\0),l). (24) 

A simple substitution 7 = (a — /3)^/r] brings us to 

AW w (a,P, VV) = f d 2 7P(a- 7 /^)AWy| 7 |,|0),l). 

(25) 

Finally, we introduce a scaled radial marginal of the P func- 
tion, 



R a (b) = - P(a-be^/^/rj)d<p, (26) 
V Jo 

which is defined for nonnegative real arguments b. For classi- 
cal states, this quantity fulfills all requirements of a probability 
density. The objective function (|24|) is then expressed by 



POO 

AW w (a,p,r))= / R a {b)AW w (b,\0),l)db. (27) 
Jo 

Until now, maximization AW H .(a, p,rj) over all classical 
states p c i is equivalent to maximization over all probability 
densities R a {b). 

We would like to add one additional constraint on R(b). 
Let us assume that the probability for finding more than N 
photons is less than some constant, i.e. 

00 N 

£ Pn (a,t) = l-Y / Pn(a,t)<l-p I . (28) 



n=N+l 



n=0 



This accounts for the fact that we did not observe more than N 
photons in an experiment, therefore the probability for doing 
so should be quite small. We can express this condition in 
terms of the P function of the input state as 



N 



d 2 pP{P) £ 



\(3-a\ 2n r] n 



-|/3-a| 2 77 



> Pr- (29) 



ri=0 



By applying the same considerations as before, and introduc- 
ing the abbreviation 



/v 



(30) 



n=0 



we can express this inequality in terms of the density R a (b) 
and obtain 



(22) 



R a {b)G{b)db > p T 



(31) 
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In conclusion, we want to optimize the objective function Then, the set of constraints can be formulated in the form 



Eq. (27 i with respect to the probability density R(b), satisfy- 
ing th 



le natural constraints 



R a (b)db=l, R a (b)>0, 



(32) 



as well as our additional requirement Eq. (31 1. This is a linear 
optimization problem on the convex set of admissible R a {b). 
Since the objective function as well as the constraints do not 
depend on a, the solution does not as well, and we will omit 
the index a from now on. 

Before we continue with the solution of this problem, let us 
comment on two properties of the function G(b). On the one 
hand, it can be easily seen that G(b) is strictly increasing with 
N for all fixed b > 0. On the other hand, it is decreasing with 
increasing b for all N > 0. The latter fact becomes clear when 
one calculates the derivative, 



G\b) = -2 



(33) 



which is strictly negative. 



2. Structure of the solution 

It is well-known that the maximum of a linear optimiza- 
tion problem is achieved at vertices of the admissible set of 
R{b) [32]. For this purpose, we have to examine the struc- 



ture of the set of functions satisfying the constraints ( 3 1 1 and 
( |32] >. As a first step, we analyze a discrete analogue of this set. 
Afterwards, we will apply the results to the continuous case. 



For discretisation, we replace the integrals in ( 3 1 1 and ( 32 1 
by Riemannian sums with K + 1 interpolation points, i.e. 

/ /(6)d6^^/(A&fc)A&. (34) 
Jo k=0 

By introducing discrete variables 

r k = R{Abk)Ab, g k = G(Abk) (35) 
with indices < fc < K, the constraints are replaced by 



A" 



^2 rfe 5fc > P 

fc=0 
K 



1. 



r k > 0. 



(36) 
(37) 



fe=0 



The examination of such convex sets of {r k ) k=0 is a standard 
problem. The vertices can be determined as follows: First, we 



rewrite the inequality ( 36 1 as an equality with an additional 



constrained slack variable y, 

K 

^2r k g k - y= Pl , 

k=Q 



y>0- 



(38) 



1 

9o 



. 1 

■ 9K 

r k > 0, 







r 

y>0 



(39) 
(40) 



To find a vertex, we set an arbitrary selection of K compo- 
nents of the vector (r, y) T to zero, and solve the equation for 
the remaining two components. If the resulting solution has 
only nonnegative entries, it is a vertex of the polyhedron of 
feasible solutions of the optimization problem. We have to 
distinguish two cases: 

1. We set K elements of the vector f to zero. Therefore, 
we are left with some element r k and the slack variable 
y, and consider the equation 



1 

9k -1 



The solution of this equations reads as 

rk = 1) Vk=gk-Pr- 



(41) 



(42) 



Nonnegativity of y has to be taken into consideration, 
therefore g k = G(Abk) should be larger than p r . Since 
this sequence of numbers is strictly decreasing with in- 
creasing fc, the index k has to be sufficiently small. 

2. We set y = and K — 1 components of the vector r to 
zero. Hence, we are left with two components r k and 
ri, k < I. This leads to 



1 1 

9k 9i 



and its solution is given by 



Tk = 



Pr ~ 91 

9k - 9i 



ri = 



gk - Pr 

9k~ 9i 



(43) 



(44) 



Since the sequence g k = G(Ab k) is decreasing, the de- 
nominator is strictly positive for pairs k < I. Therefore, 
we have to require that / is sufficiently large such that 
p T > gu and k is sufficiently small such that g k > p r . 

In conclusion, we have found that only two different types 
of vertices exist: The one only have a single element of the 
sequence r k being different from zero, provided that k is suf- 
ficiently small. The others have two elements r k , ri being dif- 
ferent from zero. We emphasize that this is independent of the 
number of sampling points K. Therefore, we may improve 
the accuracy of the approximation by an Riemannian sum by 
increasing the number of sampling points and decreasing their 
distance A&. 



3. Solution of the continuous optimization problem 

Knowing the vertices of the discretized admissible set, let 
us return to the original problem. We do this by taking the 
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limit Ab —> 0, K —> oo. The approximation of an integral by 
a Riemannian sum effectively replaces the function R(b) by a 
step function 



K 



RAb{b) = ^ r k rect Ab (b - Abk), 



(45) 



fc=0 



with rectAb(a;) = 1/A6 for all < x < Ab and rectAi>(aO 
elsewhere. All rectangles have unit area, satisfying 



/oo 
rect A h(a;) = 1. 
-oo 



(46) 



for all positive Ab. The vertices with a single nonzero com- 
ponent then read as 



R^ Ab (b)=Tect Ab (b-Abk). 



(47) 



In the limit of Ab — > 0, the function rectAb(^) converges to 
the ^-distribution, centered at some point b\. Therefore, the 
vertex takes the form 



R^(b) = 6(b-h). 



(48) 



The constraint g k > Pr has to be replaced by G(b\) > p r . 
Analogously, the vertices with two nonzero components , r; 
read as 

R v.Ab( b ) = ^^rect Ab (6 - Abk) 



9k - 9i 

9k~ Pi 

9k - 9i 



rect A fe(&- Abl). (49) 



In the continuous limit, we have to replace the rectangular 
functions by <5-functions and the coefficients — G(Abk) 
by G(b\ 2 ). Hence, these vertices are given by 



„(2) ,,s Pr-G(b 2 ) - , 
R br,b 2 \ b > = TvTl 7T7T^ 5 ( b ~ b ^ 



G{b x ) - G(b 2 ) 
G{b 1 )- Pl 



G^) - G{b 2 ) 
where the parameter b, , b 2 satisfy 

G(b 2 ) < Pr < G(bi) 



5(b-b 2 ), (50) 



(51) 



In conclusion, we have demonstrated that the objective 
function (27 1 is maximized for some function R(b) having the 



form (48 1 or ( |50| l. This solution only depends on one real pa- 
rameter 61 or two real parameters b, , b 2 respectively. There- 
fore, in order to find the systematic error, one has to find the 
global maxima of the two objective functions 



AWi 1 \b 1 )=t 2 AW w (b 1 ,\0),l) 



(52) 



AWS 2 \h,b 2 )=t 



+ 



Pr ~ G(b 2 ) 

G(h) - G(b 2 ) 
G(h)-p 



AW w (6i,|0),l) 



G(h) - G(b 2 ) 



AW w (b 2 ,\0),l) 



, (53) 



with G(b 2 ) < p r < G(bi). The largest solution gives the 
systematic error AW W . C \.. 
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